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The decay modes J/ip — ► ^7r ,^r] and 77/ are analyzed using a data sample of 
58 million J/ip decays collected with the BESII detector at BEPC. The branching 



fractions are determined to be: Br(J/^ -► 771-°) = (3.13±g;f|) x 10~ 5 , Br(J/ip -»■ 
71?) = (11.23 ± 0.89) x 10~ 4 , and Br{J/^ -► 77/') = (5.55 ± 0.44) x 10~ 3 , where the 
errors are combined statistical and systematic errors. The ratio of partial widths 
T(J/ip — > jr)')/T(J/ijj — ► 777) is measured to be 4.94 ± 0.40, and the singlet-octet 
pseudoscalar mixing angle of 77 — 77' system is determined to be Op = (— 22.08±0.81)°. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Cs 



I. INTRODUCTION 

In flavor-<S'C/(3), the 7r°, 77 and 77' mesons belong to the same pseudoscalar nonet. 
The physical states 77 and 7/ are related to the S77/(3)-octet state r/ 8 and the SUf(3)- 
singlet state rji, via the usual mixing formulae: 

i] = T] 8 cos 6p — r/i sin Op, 

rf = i] 8 sin Op + 771 cos Op, 

where Op is the pseudoscalar mixing angle [1, 2]. The conventional estimate of rj — 7/ 



mixing uses the quadratic mass matrix 
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where M| 8 = ^(Am 2 K — m 2 ) is given by the Gell-Mann-Okubo mass formula. Diag- 
onalization of this matrix gives 
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With a linear mass matrix and the linear Gell-Mann-Okubo mass formula M 88 = 
^(4m K — m n ), Op is computed to be about —24° [2]. 

The mixing angle has been measured experimentally in different ways, and the 
value is around —20° [2]. One of these measurements is based on J/ip radiative 
decays. In the limit where the OZI rule and SUf(3) symmetry are exact, one gets [3] 

R = WPI^^ = (^?- cot 2 Op, 
r(j/^ 777) y Pv J 
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where p v and p v > are the momenta of 77 and rj' in the J/ip Center of Mass System 



The first-order perturbation theory [4, 5] expression for the partial width 
T(J/ip — > 7 + pseudoscalar) is 



Here Rj/$(0) and Rp(0) are the wave functions at the origin of the J/ip and the 
pseudoscalar with mass M P) and Q c is the charge of the charmed quark. The 
pseudoscalar helicity amplitude H p (x) depends on x = 1 — {jff^) 2 ] numerically 
x\H p (x)\ 55 for Mp = m^/. Rj/^(0) and -Rp(O) can be determined from the 
J/ip — > e + e~ and P — > 77 partial decay widths, respectively. Using the lowest-order 
QCD formula for a s , the J/^ — > 777' decay width is calculated to be 213 eV, which 
is in agreement with the experimentally measured value. The value of F(J /ip — > 777) 
determined from the same formula disagrees with measurements. Some models that 
assign a small admixture of 77 and rj' to other states have been proposed to explain 
the large value of the ratio R = T(J/vp — > 777') /T( J/-?/; — > 777). For example, Ref. [7], 
which assigns small cc contribution from r] c in the 77 and 77' wave functions, predicts 
R = 3.9; Ref. [8] gives a value of R — 5.1 by considering some admixture of the 
t(1440) to the 77 and 77'. A precision measurement of R could distinguish between 
these mixing models, as well as provide a determination of the mixing angle 9 P . 
Experimental measurements of Br(J/ip — > 777) and Br(J/ip — > 777') were reported 
by the DESY-Heidelberg group [9], the Crystal Ball [10], Marklll [11] and DM2 [12]. 

The decay J/ip — > 771- is suppressed because the photon can only be radiated 
from the final state quarks. This branching fraction was measured by DASP [13] 
and Crystal Ball [10]; the average of the measurements, (3.9 ±1.3) x 10 -5 [6], is in 
agreement with the VMD prediction 3.3 x 10~ 5 [14]. In contrast, the QCD multipole 
expansion theory [15] predicts a value of 1 x 10~ 6 . 

In this paper, J ftp — > 771- is studied using 7r° — > 77 decay, J/ip — > 777 is measured 
using 77 — > 77 and 77 — > 7r°7r + 7r~ with tt° — > 77, and J/t/> — > 777' is studied using 77' — > 
77, 77' — * 77r + 7r~ and 77' — > 777r + 7r~ with 77 — > 77. The analyses use a data sample 
that contains 58 x 10 6 J/ip decays collected with the updated BEijing Spectrometer 



(CMS). 
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(BESII) operating at the Beijing Electron Positron Collider (BEPC). 
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II. BES DETECTOR AND MONTE CARLO SIMULATION 

BESII is a large solid-angle magnetic spectrometer that is described in detail in 
Ref. [16]. The momentum of charged particles is measured in a 40-layer cylindrical 
main drift chamber (MDC) with a momentum resolution of <r p /p=1.78%\/l + p 2 (p 
in GeV/c). Particle identification is accomplished using specific ionization (dE/dx) 
measurements in the drift chamber and time-of-flight (TOF) information from a 
barrel-like array of 48 scintillation counters. The dE/dx resolution is CdE/dx — 
8.0%; the TOF resolution for Bhabha events is (Jtof = 180 ps. Radially outside 
of the time-of-flight counters is a 12-radiation-length barrel shower counter (BSC) 
comprised of gas tubes interleaved with lead sheets. The BSC measures the energy 
and direction of photons with resolutions of a E /E ~ 21%\/E (E in GeV), a<f, = 7.9 
mrad, and o z = 2.3 cm. The iron flux-return of the magnet is instrumented with 
three double layers of proportional counters that are used to identify muons. 

A GEANT3-based Monte Carlo simulation package [17], which simulates the de- 
tector response including interactions of secondary particles in the detector material, 
is used to determine detection efficiencies and mass resolutions, optimize selection 
criteria, and estimate backgrounds. Reasonable agreement between data and MC 
simulation is observed for various calibration channels, including e + e~ — > (j)e + e~, 
e + e~ — > (7)/i + /i~, J/ip — > pp, J/ip — > pn and ip(2S) — > n+Ti^J/ip, J/ip — > l + l~ . 

III. DATA ANALYSIS 

A. J/V> - * 777 

In the J/ip — > 777 decay mode, there are no charged tracks in the final states. 
Each candidate event is required to have three and only three photon candidates; 
the MC indicates that the number of these decays that produce final states with 
more than three photon candidates is negligible. A photon candidate is defined as 
a cluster in the BSC with an energy deposit of more than 50 MeV, and with an 



angle between the development direction of the cluster and the direction from the 
interaction point to the first hit layer of the BSC that is less than 20°. If two clusters 
have an opening angle that is less than 10° or have an invariant mass that is less 
than 50 MeV/c 2 , the lower energy cluster is regarded as a remnant from the other 
and not a separate photon candidate. A kinematic fit that conserves energy and 
momentum is applied to the three photon candidates, and x 2 < 20 is required. We 
also require |cos#u| < 0.8 and 6 m i n > 6°, where 9 V is the polar angle of a decay 
photon in the pseudoscalar's CMS (shown in Fig. la), and Q m in is the minimum 
angle between any two of the three photon candidates (shown in Fig. lb). This 
rejects background from the continuum e + e~ — > 77(7) process. 
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FIG. 1: Distribution of (a) cos 6^ and (b) 6 m i n . The open histograms are J/ip data, the 
shaded histograms are background from e + e~ — ► 77(7), and the dashed lines are simulated 
J/ip — ► 77T° — ► 777 events (not normalized). 



1. J/lp — > 77T°, 7T° — > 77 

Figure 2 shows the invariant mass distribution in the 7r° mass region of the two 
photon candidates that have the smallest opening angle. A peak at the n° mass is 
evident. 

From MC studies, background channels that produce a peak in the 7r° signal 
region come mainly from channels with 57 final states, such as J/ip — » 77r°7r°, via the 
/ 2 (1270), / (2100) etc. (J ftp — * 47s violates C-parity). These background sources 
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FIG. 2: Invariant mass distribution of the 77 with smallest opening angle for J/ip — ► 777 
candidate events. The solid squares with error bars are data, the histogram is the best fit 
described in the text, and the dashed line is the background. 

are studied using events where the number of photon candidates in the event is four. 
Four photon events are selected and subjected to a four-constraint kinematic fit to 
J/ip 777, using any three of the four photons; the three-photon combination 
with the smallest x 2 is selected for the background study. Figure 3 shows the 
invariant mass distribution for the two photons with the smallest opening angle 
from four-photon events. A peak is observed in the 7r° mass region that agrees with 
expectations from MC simulations that include all known modes that produce 57 
final states. However, since the known background channels do not account for the 
level of the observed background in the data sample, a scale factor is introduced to 
scale the MC background predictions for fits to the distribution in Fig. 2. The scale 
factor depends strongly on which intermediate states are considered for J/ip — > 57 
decays; the difference between the scale factors determined from different channels 
is treated as a systematic uncertainty of the background subtraction. 

Figure 2 is fit with a MC-simulated J/ip — > 771- histogram for the signal, a MC- 
simulated J ftp — > 57 background shape, as well as a shape of MC simulated phase 
space for other sources of backgrounds. The number of 771° events determined from 
the fit is 586 ±51. The MC-determined detection efficiency for J/ip — > 771- , ir° — > 77 
is e = (32.80 ± 0.21)%, where the error comes from the limited statistics of the MC 
sample. 
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FIG. 3: The invariant mass distribution of 77 pairs with the smallest opening angle in 
J/ip — ► 777 events selected from the four photon event sample. 

2. J/ip — > 777,77 — > 77 

Figure 4 shows the invariant mass distribution of the two photon candidates with 
the smallest opening angle in the r\ mass region, where an r\ peak is evident. 
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FIG. 4: Invariant mass distribution of the 77 with the smallest opening angle of J/ijj — > 
777 candidates. Solid squares with error bars are data, the histogram is the fit result, and 
the dashed line is the background. 

The 77 invariant mass distribution of Fig. 4 is fit with a histogram from MC- 
simulated J/ip — > 777, 77 — > 77 events and a second order Legendre polynomial back- 
ground function. The fit yields a signal of 9096±133 r/s. The MC-determined 
detection efficiency is e = (36.33 ± 0.22)%. 
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3. J/tp — > 7^') 77' — > 77 

Since the momentum of the rj' is lower than that of the 7r° and 77 in J/ift radiative 
decays, the angle between the two rj' decay photons is not small enough to be 
useful for distinguishing them from the radiative photon. For this channel, the mass 
distribution of the three 77 combinations for each event are plotted in Fig. 5, where 
an rj' signal is evident above a smooth background due to wrong 77 combinations 
plus other background sources. 
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FIG. 5: The 77 invariant mass distribution for J/ip — > 777 candidate events (three entries 
per event). The solid squares with error bars indicate data, the histogram is the fit result, 
and the dashed line is the non-combinatorial background. 

A fit to the data points, with the MC simulated mass distribution for the 
J/ijj — > 777', r( — > 77 decay including combinatorial background for the signal and 
a second order Legendre polynomial for background between 0.8 and 1.2 GeV/c 2 , 
yields 2982±101 entries. Since all the 77 combinations are plotted in the M 77 distri- 
bution, the combinatorial background is included in the entries for both data and MC 
simulation. The efficiency for signal J ftp — > 777', rj — > 77 entries is (40.30 ± 0.22)%. 
The combinatorial background is about 20% for both data and MC simulation, and 
they cancel out when the N obs is divided by the efficiency 40.30% in the branching 
fraction calculation. 
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B. J ftp — > 7777T + 7T 

In this final states, there are two charged particles 7r + and n~ and three photons. 
Candidate events are required to satisfy the following common selection criteria: 

1. Two good charged tracks with net charge zero. Each track must have a good 
helix fit, a transverse momentum larger than 60 MeV/c, and |cos#| < 0.8, 
where 8 is the polar angle of the track, and must originate from the interaction 
region. 

2. At least one charged track is identified as a n, satisfying Xpioi 71 ) < Xpid(K) 
and Xpid^) < Xp/d(p)> where X 2 pid = X 2 dE /dx + Xtof is determined using 
both dE/dx and TOF information. 

3. At least three photon candidates are required. The photon identification is 
similar to that used in the J/ift — > 777 analysis, except that the angle between 
a cluster and any other cluster must be greater than 18°, and the angle between 
the cluster and any charged track must be greater than 8°. These differences 
reflect different sources of fake photons. 

4. A four-constraint kinematic fit is applied to all three-photon combinations 
plus the two charged tracks assuming J/ip — > 7777r + 7r~. The three-photon 
combination with the smallest x 2 is selected, and the \ 2 of the kinematic fit 
is required to be less than 20. 

The events that survive these selection criteria with an invariant mass in the 
range M^+j- < 1.2GeV/c 2 are assumed to come from either 77 or rj' decays, and the 
other photon is considered to be the radiative photon. Figure 6 shows a scatterplot 
of M 77 versus M 117T + n - for the selected events. Clear rj and rj' signals corresponding 
to r] — > 7r°7r + 7r~, 7T° — > 77, and rj 1 — > r)7r + ir~ , rj — > 77 are observed. 

1. J/tp — ► 77/, 77 — > 7r°7r + 7r~ 

After the requirement that the 77 invariant mass is in the 7r° mass region (M 77 e 
[0.088, 0.182] GeV/c 2 , ±3<r), a clear 77 signal is evident in the 777r + 7r~ invariant mass 
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FIG. 6: Scatterplot of M 77 versus M yyn + n - for the J/ip — > 7777r + 7r candidates, 
distribution shown in Fig. 7. 
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FIG. 7: The ^^tt + tt~ invariant mass distribution for J/ip — > 7777r + 7r~ candidates that 
satisfy the requirement M 77 G [0.088,0.182] GeV/c 2 . The solid squares with error bars 
indicate the data, the histogram is the fit result, and the dashed line is the background. 

The simulated M 117T + W - mass distribution from the signal MC and a second-order 
Legendre polynomial are used to fit the 777r + 7r~ invariant mass distribution. The 
fit gives 1885 =b 58 77 events. The MC-determined detection efficiency for J/ip — > 
-fT],V^ vr°7r+7r-, and tt° -> 77 is e = (12.25 ± 0.15)%. 



2. J/ip — ► 7V) V — * W +7r 

The 777r + 7r~ invariant mass distribution for events with 77 mass within 3a of 
the 77 mass (M 77 e [0.484,0.612] GeV/c 2 ), is shown in Fig. 8. 
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FIG. 8: The 777r + 7r~ invariant mass distribution for events with 77 mass in the 77 mass 
region (M 77 G [0.484,0.612] GeV/c 2 ). The solid squares with error bars indicate data, the 
histogram is the fit result, and the dashed line is the background. 

A similar fit as for 77 — > 7r°7r + 7r~ yields 8572 ± 131 7/ events; the MC-determined 
detection efficiency for J/t/j — > 777', 77' — > r]iT + ir~, and 77 — > 77 is £ = (16.10 ± 0. 12)%. 



C. J/tp —>■ 777T + 7T 

J/V — > 777' is also studied using the 77' — > 77r + 7r~ decay channel. For this study, 
the 7r ± and photon selection requirements are the same as used for the J/ijj — > 
7777r + 7T~ final state, and the event selection is similar, except that here at least two 
photons are required in the event. The photons and charged tracks are kinematically 
fitted to J/ip — > 777r + 7r~ assuming four-momentum conservation, and x 2 < 20 is 
required. When there are more than two photons, the kinematic fit is repeated 
using all possible photon combinations, and the one with the smallest x 2 is kept. 
The photon with the higher energy is considered to be the radiative photon from 
the J /ip decay. Figure 9 shows the invariant mass distribution of 77r + 7r~ for the 
candidate events where an 77' signal is evident. 

Figure 9 shows the result of a fit to the 77r + 7r~ invariant mass distribution that 
follows a similar procedure as that for the fit to the 777r + 7r~ distribution of the 
previous section. The fit yields 23243 ± 229 77' signal events. The MC-determined 
detection efficiency for J/ip — > 777', 77' — > 77r + 7r~ is e = (25.02 ± 0.10)%. 
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FIG. 9: The 77r + 7r - invariant mass distribution for selected J ftp — > 777r + 7r~ events. The 
solid squares with error bars indicate data, the histogram is the fit result, and the dashed 
line is the background. 

IV. SYSTEMATIC ERRORS 



Systematic errors in the branching fraction measurements mainly originate from 
photon identification (ID), MDC tracking efficiency, particle ID, kinematic fitting, 
mass resolution, 7r° reconstruction, and parameterizations of background shapes. 



A. Photon identification 

The efficiency for photon ID is discussed in Ref. [18]. It is found that the rela- 
tive efficiency difference between data and MC simulation for high energy photon 
detection is about 0.8% per photon, while for low energy photons, the difference is 
around 2% per photon. Since the energy of the radiative photon in J/ip — > 77r°,7?7, 
and 777' is high, and the energies of the photons from pseudoscalar particle decays 
are low, the total systematic error due to photon ID is taken as (0.8 + 2. On) %, where 
n is the number of photons from the pseudoscalar particle decay. 

B. MDC tracking 

The MDC tracking efficiency is studied in Ref. [17]. It is found that there is a 
2.0% relative difference per track between data and MC simulation. For the channels 
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in this analysis that have two charged tracks, a 4% systematic error on the MDC 
tracking efficiency is assigned. 

C. Particle ID 

A clean charged n sample obtained from J ftp — > pn without the use of particle 
ID is used to study data-MC differences between particle ID efficiencies for different 
momentum ranges. Since only one of the two charged tracks is required to be 
identified as a pion, the MC simulates data rather well; it is found that the MC 
simulation agrees with data within 0.2% for both J/ip — > 7777r + 7r~ and J/xjj — > 
777r + 7r~ modes. 

D. Kinematic fit 

Samples of J/ip — > pn and e + e~ — > 77 events selected without using kinematic 
fits are used to study the systematic error associated with the four-constraint kine- 
matic fit. For the x 2 < 20 criteria, the difference of kinematic fit efficiencies between 
data and MC simulation is less than 1.2% for pir, and 2.4% for e + e" — > 77. Ex- 
trapolating these differences to the channels reported here, we conservatively assign 
a 4% systematic error to the kinematic fit efficiency. 

E. Different mass resolution between MC and DATA 

There is a slight difference of the mass resolution between MC simulation and 
data. When the histogram shape of invariant mass distribution from MC simulation 
is used to fit the invariant mass distribution of data, it introduces some systematic 
error. The high statistics decay channels J/ip — > 777,77 — > 77, J/ip — > 7V) v' ~^ 
777T+7T - and J/tp — > 777', rj' — > r yn + ir~ are used to study this source of systematic 
error. For these channels, we allow the mass resolution to vary in the fit to the 
invariant mass distributions, and we also determine the number of signal events 
by subtracting side-band-estimated backgrounds. The resulting branching fractions 
change by at most 1.6%, 0.1%, and 0.6% for J/tp — > 777,77 — > 77, J/ip — > 777', r/' — > 
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r]iT + iT~, and J ftp — > 777', 77' — > 77r + 7r _ respectively. We assign 1.6%, 0.1% and 0.6% 
as the systematic errors due to mass resolution uncertainties for the J/ip — > 777, 
J/ip — > 7777r + 7r~ and J/^ — > 777r + 7r~ decay modes, respectively. 

F. Reconstruction of 7r° 

In J/^ — > 77r°, the 7r° momentum is high and the angle between the two decay 
photons is small. As a result, it is possible for the two photons to merge into a 
single BSC cluster. According to a study reported in Ref. [19], the systematic error 
associated with 1.5 GeV 7r° reconstruction is 0.83%. The effect on low energy 7r°s 
or 77s is small enough to be neglected. 

G. Background shape 

For the J/ip — > 771° mode, the background estimate based on the four photon 
event sample has a large uncertainty. Fits using MC-determined background shapes 
from different background channels yield different numbers of signal events; the 
corresponding changes in the branching fractions range between tg 6 g 4 %. The largest 
difference is taken as the systematic error. Different order Legendre polynomials 
are used to fit the mass spectra for the other decay modes, and the differences 
between these fits and those used to get the numbers of signal events are used as 
the systematic error due to background parameterization. Different fitting ranges 
are also used in the fit, and the differences are included in the systematic error. The 
uncertainty due to the background shape and fitting range is less than 2%. 

H. Branching fractions of the secondary decays 

The branching fractions of decay from n°, rj and 77' are taken from the PDG [6]; 
the uncertainties are included in the measurement errors of the reported branching 
fractions. 
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I. The number of J ftp events 

The total number of J/ip events, determined from the 4-prong data sample, is 
(57.7 ± 2.72) x 10 6 . The 4.72% relative error is taken as a systematic error [20]. 

J. Total systematic error 

Table I summarizes the systematic errors from all sources for each mode. We 
assume all the sources are independent and add them in quadrature; the resulting 
total systematic errors are ±™i%, 8.1%, 10.6%, 9.3%, 9.5%, and 8.7% for J/ip -> 
77r° — > 777, J/ip — > 777 — > 777, J/tp — > 77/ — > 777, J/if) —> jr) —> 7777r + 7r~, 
J/^ — > 7?]' — > 7777r + 7r~, and J/^ — > 777' — > 777r + 7r~, respectively. 

TABLE I: Summary of the systematic errors (%). 



Sources 


7T° —> 77 7] 


77 rj 


' -► 77 n 


— ► 7T 7r + 7T 


rj' — > r]TT + ir 


T]' — > 77T+7T 


Photon ID 


4.8 


4.8 


4.8 


4.8 


4.8 


2.8 


Tracking 








4.0 


4.0 


4.0 


Particle ID 








0.2 


0.2 


0.2 


Kinematic fit 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Mass resolution 


1.6 


1.6 


1.6 


0.1 


0.1 


0.6 


7T° reconstruction 


0.83 












Background shape 


+16.4 
-6.8 


0.73 


1.8 


1.7 


0.5 


0.2 


Branching fraction used 


0.04 


0.66 


6.61 


1.77 


3.45 


3.39 


Number of J ftp 


4.72 


4.72 


4.72 


4.72 


4.72 


4.72 


Statistic of MC sample 


0.64 


0.61 


0.55 


1.23 


0.75 


0.40 


Total error 


+18.3 
-10.6 


8.1 


10.6 


9.3 


9.5 


8.7 
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V. RESULTS AND DISCUSSION 

The branching fractions of J/ decays are determined from the relation 

BrU/i, - 7 P) = N^VN^lJl^jn 

K IV 1 ' N J ^ ■ Br(P -^Y)-e{J/^ ^ 1 P -> 7 y) ' 

where P is either 7r°, 77, or 77', Y is the pseudoscalar decay final state, and Br(P — > F) 
is the branching fraction of the pseudoscalar decays into final state Y. The results 
of Br{J/ip -> 7P) are listed in Table II. 

The branching fractions of J/ip — > 777, J/^ — > 777' measured from different decay 
modes are consistent with each other within the statistical fluctuations and uncom- 
mon systematic errors. The measurements from the different modes are, therefore, 
combined using a standard weighted least-squares procedure taking into considera- 
tion the correlations between the measurements; the mean value and the error are 
calculated by: 



1 



Here u>ij is the element of the weighted matrix W = V~ x , where V x is the covariance 
matrix calculated according to the systematic errors listed in Table I. For J/ip — > 777, 
the correlation coefficient between 77 — > 77 and 77 — > 777r + 7r~ is p(l,2) = 0.553; 
for J/ip — > 777', the correlation coefficients between 77' — > 77, 77' — > 77r + 7r _ and 
77' -> 777T+7T- are p(l, 2) = 0.296, p(l, 3) = 0.404 and p(2, 3) = 0.703. The weighted 
averages of BESII measurements and the PDG [6] values are listed in Table II. 

Figure 10 shows the comparisons between the measurements in this paper and 
those from previous measurements [9-13]. Our measurement of Br(J/ip — > 771" ) 
agrees with those of Crystal Ball [10] and DASP [13] within the large errors of 
the previous measurements, and has much improved precision. Our measurement's 
lower central value may be because background channels that produce a peak in 
the signal region have been considered. Our measurements of Br(J/ip — > 777) and 
Br(J/ip — > 777') are higher than the PDG world average [6], and have better precision 
than the previous measurements [9-12]. 

The results listed in Table II also allow us calculate the relative branching frac- 
tions for 77 and 77' decays; considering the common errors in the measurements, one 
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TABLE II: Branching fractions of J/tp — ► 771"°, 777 and 77/'. 



Decay mode 


BESII 


BESII combined 


PDG [6] 


7 7r° 


7T° — * 77 


(3.13 ± 0.28±g:||) x lO" 5 


(3.i3±g:S|) x 10- 5 


(3.9 ± 1.3) x lO" 5 


777 


77 — ► 77 


(11.00 ±0.16 ±0.90) x 10~ 4 


(11.23 ±0.89) x 10~ 4 


(8.6 ±0.8) x 10~ 4 




77 — * 7r°7r + 7r~ 


(11.94 ± 0.37 ±1.11) x 10~ 4 








-> 77 


(6.05 ±0.21 ±0.65) x 10~ 3 






77?' 


7/ -> 7p 


(5.46 ±0.06 ±0.48) x 10~ 3 


(5.55 ±0.44) x 10~ 3 


(4.31 ±0.30) x 10~ 3 




7/ — > 7y7r + 7r _ 


(5.28 ±0.08 ±0.51) x 10~ 3 











BESII 








CBAL 








DASP 

1 






5 


10 


15 



10- 













BESII 






CBAL 






DASP 






CNTR 















BESII 


6 














MARKIII(1) 


4 






MARKIII(2) 








_ DM2(1) 


2 






DM2(2) 








CBAL 












11 13 15 17 19 
lO" 4 



2.5 



5 

lO" 3 



7.5 



10 



FIG. 10: Comparisons of Br(J/ip — ► 7vr°), Br(J/ip — > 777) and Br(J/t/j — > 777') between 
BESII and previous measurements [6]. The shaded regions are the world averages from 
the PDG [6]. 



gets 



Br{r}' — > 77) 



Br(J/ip -> 7»7'>»7' -> 77) 



Br(j]' — ► 77T+7T ) _Br( J/?/; 777', 77' 77T+7T ) 

Br(rj' — > r)ir + Tr~) Br{J ftp — > 777', 77' — * -i]tt + tt^) 

Br(rj' — > 77r+7r _ ) Br{J/il> — > 777', 77' — > 77r+7r~) 

Br(?7 — > 77) Br(J/ip — > 777, 77 — + 77) 



0.080 ±0.008, 



= 1.45 ±0.07, 



1.61 ±0.14. 



.Br (77 — > 7T°7r+7r~) Br(J/ip — > 777,77 — * 7r 7r+7r~) 

The correlation coefficients between denominator and numerator in above equations 
are 0.419, 0.859 and 0.575 respectively. The world averages [6] of the same ratios 
are 0.072 ± 0.006, 1.50 ± 0.08 and 1.75 ± 0.04 respectively. The agreement is quite 
good. 
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If both the OZI rule and the S77/(3) symmetry are exact, it is expected that [3]: 

r(j/^ T ,o ^ 
r( j/*. -. 7 ,) ^ p, 

Using Br(J/ip — > 777) and Br(J /ip — > 777') in this analysis, one obtains 

i? = 4.94 ± 0.40, 

|0 P | = (22.08 ±0.81)°, 

where the common errors have been considered in the ratio calculation. Comparing 
with the mixing models with states other than 77 and 77', the measurement of R 
agrees with the prediction of R = 5.1 [8] within one standard deviation, while it 
deviates from R = 3.9 [7] by more than 3 standard deviations. According to the 
theoretical calculation of Ref. [2], the value of 6p is negative, in which case its value 
is 6 P = (-22.08 ±0.81)°. 



VI. SUMMARY 



Using 58 million J /if) events collected by BESII, the branching fractions of J/ip 
decays into a photon and a pseudoscalar meson are measured as Br(J/ip — > 777°) = 
(3.13t8i4) x 10~ 5 , Br(J/ip -> 777) = (11.23 ± 0.89) x 10~ 4 , and Br(J/ip -> 777') = 
(5.55 ± 0.44) x 10~ 3 . The results are compared to 77 and 77' mixing models. 
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